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INTRODUCTION 

An optimized cost and performance model for a phosphoric acid fuel cell 
powerplant system has been derived and developed into a modular FORTRAN 
computer code. The mathematical model developed in this study combines 
appropriate cost, energy, and electrochemical analysis as applied to the major 
components of a phosphoric acid fuel cell system. 

In the model, the following important parameters were optimized: steam to 

methane ratio in the reformer, hydrogen utilization in the fuel cell stack, and 
the number of fuel cell plates per stack. 

The optimization technique used here was by means of COMPUTE code which 
was developed for soliving nonlinear programming problems. 


I. SYSTEM DESCRIPTION 

As shown in Figure 1, methane which is circulated by compressor (C) is 
preheated by heat exchanger E-1 prior to mixing it with the super heated steam 
which receives its beat by passing through heat exchanger E-9. Before entering 
the reformer, the methane steam mixture is heated via heat exchangers E-2 and 
E-3. Inside the reformer, methane is catalytically reformed by reaction with 
excess steam to produce carbon monoxide, carbon dioxide, and the desired pro- 
duct, hydrogen. The effluent from the reformer is cooled by flowing through 
heat exchanger E-2 before it enters the high temperature shift converter S-1. 
The function of the high temperature shift converter is to increase the hydro- 
gen concentration and to reduce the carbon monoxide concentration of the 
reformer gas effluent. The temperature of the effluent from the shift conver- 
ter S-1 is then reduced by passing through heat excangers E-1, E-9 and E-6 
before entering the low temperature shift converter S-2. The low temperature 
shift converter further increases the hydrogen concentration by promoting the 
shift reaction at a lower operating temperature. The effluent from the low 
temperature shift converter then enters the fuel cell containing H2, CO, CH4, 
C02 and H20. The fuel cell converts inputs of hydrogen and oxygen to DC power, 
water and heat. Oxygen is delivered to the fuel cell by air compressor A, 
which also provides air to the reformer burner. The spent fuel from the fuel 
cell anode goes to the burner after mixing with air supplied by compressor A. 

Before entering the burner, the mixture is preheated by the burner 
effluent via heat exchanger E-4. The spent fuel is then burned with whatever 
additional methane is needed to provide the thermal energy necessary for the 
reformer reaction. 
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Figure 1 Flow dlagru of CSU designed PArc qrateai 














Heat generated in the fuel cell is removed by heat exchangers E~7 and E-10. 
Heat from heat exchanger E-7 can then be utilized in industrial heat processing 
or space heating and cooling, while exchanger E-10 is used to preheat the water 
supplied by liquid separator Q to provide the necessary steam needed for the 
reforming process. The effluents from the burner and fuel cell cathode will 
have their water removed and separated by condenser E-5 and liquid separator Q 
before allowing them to be exhausted to the atmosphere. 
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II. OPTIMIZATION MATHEMATICAL MODEL 

This chapter describes a practical procedure for the optimization in 
designing a processing system like the 100 kW phosphoric acid fuel cell (PAFC) 
powerplant shown in Figure 1. The design criteria are the investment cost and 
the operation cost, both to be minimized. 

Computation time increase rapidly with the number of variables to be 
optimized. The optimization procedure described here will, for brevity, be 
limited to three of the most important variables. Therefore, we will describe 
the procedure for the determination of the optimum design of hydrogen utiliza- 
tion in the stack, the steam to carbon ratio in the reformer, and the number 
of cell plates per stack. The derivation of the optimization model was based 
on the following assumptions. 

1. All reactants and products are assumed to be ideal gases. 

2. All heat exchangers in the PAFC system are assumed to be of the counter 
flow type. 

3. In the reformer and shift converters, all processes are assumed to be 
isothermal. In addition, only the demethanation and water shift reactions 
(Ref. 1) occur in the reformer, and the latter reaction reaches the equilibrium 
condition. 

4. Outlet temperature of the fuel cell stack is equal to the stack 
operating temperature. 
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5, The combustion process in the burner is assumed to go to completion. 


To carry out the required optimization, performance and cost analyses of the 
PAFC system will be essential. Energy, mass, and electrochemical analyses in 
the reformer, the heat exchangers, the shift converters, and the fuel cell 
stack were combined to develop a mathematical model and the associated computer 
program for PAFC powerplant (Ref. 1). 

In the following cost analysis and formulation of the optimization problem 
will be presented. 

2.1 Cost Analysis 

The cost, which is to be minimized,, is the sum of the equipment capital 
costs. For the sake of simplicity, it will be assumed that the cost of piping 
and auxiliary equipment (like pump, compressors, and separator) are negligible. 

To minimize equipment cost, expressions for the installed cost for each 
primary equipment are required. Such expressions are readily obtained by 
using the exponent rule for cost estimation (Ref. 2). The rule expresses the 
investment of each unit equipment to be a function of paraments that can be 
evaluated in terms of computed quantities and decision variables. Based on 
mid-1981 money, the following are installed cost expressions for the reformer: 
Cr, high and low temperature shift converters: Chs and Cls, respectively, fuel 
cell module; Cf, and beat exchangers: Ce. 



Cr 

X 309.82 

(1) 

Chs 

X 162.65 {H7+H9/®*®^ 

(2) 

CIS 

X 230.34 (H9 + H13)®‘®^ 

(3) 

Cf 

« 8.5136 (N)(Nst) 

(4) 

Ce 

» 1046.23 (Al^*^^^^ + + ^30.6924 + ^^0.6924 

+ /\50*6924 ^ ^g0.6924j 

(5) 


where W : amount of catalyst in the reformer, lb 

Hi : flow rate of steam i, Ib-mole/hr 
N ; number of fuel cell plates in one stack 
Nst: number of stacks in the system 

Aj : transfer area of heat exchanger j, m 

The cost of the power inverter can be obtained from interpolating the following 
cost table (Ref. 3), 

Voltage-V: 50 164 203 248 304 366 433 528 657 

Cost-?/V: 200 160 150 140 130 120 110 100 90 

2.2 Formulation of the Optimization Problem 

In designing the PARC powerplant, a number of constraints were imposed on 
its construction and operation. The following is description of these limita- 
tions: 


1. Since carbon monoxide is considered as poison for the fuel cell cata- 
lyst, a limit of 1 percent volume in the fuel cell feed was placed. Thus, the 
first optimization contraint is: 
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( 6 ) 


gl » Yco < 0.01 

where Yco is the mean mole fraction of carbon monoxide. 

2. Because platinum and carbon components corrode at cell potentials 
above 0.8 V, and since lower potential means lower cell efficiency, the cell 
operating voltage is to be limited to the range of 0.5 to 0.75V. This voltage 
range provides the following constraints: 


g2 = V < 0.75 

(7) 

g3 = V < -0.5 

(8) 


3. For small fuel cell module consideration, it is assumed that the 
maximum output voltage per stack not to exceed 300 V. This provides the 
fourth constraint: 

g4 = Vs < 300 (9) 

where Vs is the output voltage per stack. 

4. The system is designed for 100 kW electric output. This provides a 
fifth constraint that can be expressed as, 

g5 = (1.0148 - [1.0298-0.001689 (5.0976-DC/lOOO)]^^^} 

/O. 0008444 < 100 (10) 

where DC is the DC electric power output. 

5. Since the electrolyte, phosphoric acid, evaporates significantly at 
temperature above 493K, and because lower operating temperature results in 
lower cell efficiency, the inlet temperature of the fuel entering from the 
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anode side and the process air entering from the cathode side is to be limited 
to the range of 398K to 477K, Thia provides the following constraints; 


G6 * T12 < 477 

( 11 ) 

G7 « T12 < -398 

(12) 

68 « T31 < 477 

(13) 

G8 « T31 < -398 

'.mm 

(14) 


where subscripts 12 and 13 refer to the stream numbers 12 and 13 in Figure 1. 

6. To provide high conversion in fuel processing subsystem, it is assumed 
that the concentration of hydrogen entering the cell anode is greater than 60 
percent. This provides the following constraint, 

GIO « -H13/F13 < -0.6 (15) 

where, H and F refer to amount of hydrogen and fuel, respectively, and sub- 
script 13 refer to the steam number 13 in Figure 1. 

With the formulation of the constraints imposed on the system design being 
completed, the optimization problem can be stated as 

Minimum ^ Ci = Cr Chs + Cls + Cf + Ce + Cv (16) 

Subject to; gl, g2, ..., glO (Equations 6 to 15) (17) 

this is a constrained minimizing problem with three decision variables and ten 
constraints. 

2.3 Optimization Algorithm 

The constrained optimization problem (equations 16 and 17) was solved by 
using mixed penalty function together with unconstrained optimization tech- 
nique for extremization. This nonlinear constrained problem, 
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( 18 ) 


Mi„. |S . 2^ Cl * K lr.(gp 

Where, K Is the penalty coefficient and ^ is a new objective function. 

A Fortran computer program had been developed to solve the optimization 
problem (Ref, 4) called COMPUTE, which will be described in the next chapter. 
Values of (ji were minimized for a sequence of monotonically increasing K > 0. 
Then, by means of Hooke and Jeeve pattern search method (Ref. 5), the optimal 
point can be determined from an assigned input for different values of the 
penalty function. 
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III. OPTIMIZATION COMPUTER MODEL 

The MAIN program calls thret subroutines which are PATTRN, OBJECT, and 
COEFF, Subroutine OBJECT calls four subroutines which are DMIXl, SCALE, RLIN, 
and HEXCl, where DMIXl and HEXCl are the same as those used in the performance 
model and RLIN is also used in the cost model. Subroutine CONVl is called by 
Subroutine BMIXl for solving a nonlinear algebraic equation. The overall 
structure of the program for evaluating some optimal design values of PAFC 
powerplant is illustrated in Figure 2. 


2.1 General Optimization Program - COMPUTE 

I ami I III, I ■ ■ ii ' .aii.iiw " ■■ « ■■■« 

The COMPUTE code (Ref. 5) was used as the general optimization program. 
COMPUTE solves the constrained optimization problems using mixed penalty 
function* together with unconstrained optimization techniques for extremiza- 
tion. There are three optimal search techniques which can solve our 
optimization problems with the mixed penalty function. These are: Hooke and 
Jeeves pattern search (Rof. 4), Fletcher and Reeves conjugate-gradient search 
(Ref. 6), and Fletcher and Powell's variable-metric search (Ref. 6). The 
Hooke and Jeeves pattern search method was selected. It performs one search 
from an input starting point. To save computation time, we also fixed the 
initial penalty level instead of using Subroutine COEFF to find one. 


(r 
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*By mixed penalty functions, we mean if the first constraints are inequal i- 
ties, and constraints (2+1) tom are equalities, our problem becomes: 


The function 0 (21, K) is then minimized for a sequence of monotonically 
increasing K > 0. 
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Pigure 2 Overall Structure of Optimization Computer Ifodsl 






There are one Main program and three Subroutines adopted from COMPUTE code 
to solve this optimization mathematical model (Equations 16 and 17). These 
will be discussed in the following sections. 

MAIN Program 

The MAIN program reads the input data and checks whether the initial 
guessed data in the limits of the assigned bounds, then controls the procedures 
for the convergence of the penalty function. The nonlir,^ar constrained problem 
is first changed to the unconstrained problem by setting up the penalty 
function. Then, by means of pattern search technique, one can find the optimal 
point at different penalty levels. For a series of increasing values of 
penalty coefficient (or penalty level), the sequence of resulting solutions 
will more and more closely satisfy the constraints. When all the constraints 
are satisfied within some preassigned tolerance, the procedure terminates. 

The nomenclature of MAIN program is shown in Table 1, while the flow chart 
of this program is shown in Figure 3. 

Subroutine OBJECT 

This subroutine computes the modified objective function. Constraints 
and objective function defined beforehand are calculated in this subroutine. 

The trial and error procedure for simulating the continuous PAFC system is 
the same as those in the performance model, and will not repeat here. The 
cost estimated in the objective function is shown in Equation 16, and the 
constraints are Equations 6 to 15. 
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TABLE 1 

NOf€NCLATURE IN MAIN PROGRAM AND SUBROUTINE OBJECT 

— - ** - ■ I I i I -^11 I — ■ ■ . M — II. II II I 

N: number of independent variables 

JEQ: number of functional equality constraints 

JINEQ: number of functional inequality constraints 

M: number of functional constraints 

IPMAX; maximum number of penalty levels 

lOBJ: function evaluation counter 

IP: penalty level counter 

IR: random check improvement counter 

IT: search iteration counter 

ITMAX: maximum number of iterations per gradient search 

ITOBJ: total number of iterations per gradient search 

ITOT: total number of search iterations 

\ NTEST: number of randomly generated exploratory points in the pattern search 

ALPHA: scaling factor between SIGMA2 and SIGMAl 

ARAND: initial parameter for the random number generator 

CFAC: proportionally constant used in calculating PNLTY 

CONST: penalty multiplier 

; CONVG: penalty level convergence criterion 

•I 

DELTA: finite-difference increment 

’ DX: length of a gradient move 

EFAC: search increment reduction factor 

, EPSLN: gradient vector convergence factor 

GRADX: gradient vector convergence criterion 

! 
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TABLE 1 (cont'd) 

NOMENCLATURE IN MAIN PROGRAM AND SUBROUTINE OBJECT 


THETA: maximum time allowed to evaluate the objective function when 

calculating the gradient vector by central differences 

GB(I): unsealed value of constraint I 

G(I): scaled value of constraint I 

SC(I): scaling factor for constraint I automatically determined at the 

beginning of each penalty level: G(I )=:GB(I)’^SC(I) 

PD(I): sum of the magnitudes of D(GB{I))/(D(X(KFR)) , KFR=:1, 2, N 

SIGMAO: sum of the squares of the constraints which are violated (excluding 

bounds on the independent variables) 

SIGMAl: sum of the squares of the scaled constraints which are violated 
(excluding the independent variable bounds) 

SIGMA2: penalty function composed of the upper and lower bounds on the 

independent variables 

SIGMA: sum of SIGMAO and a modified SIGMA2 obtained by dividing SIGMA2 by 

the penalty coefficeint squared 

TMTOT: total search time 

Program Options 

I0UT=0 output follows final optimization 

I0UT=1 output follows optimization for each starting point 

I0UT=2 output follows optimization for each penalty level 

I0UT=3 output follows each major search or gradient move 

KB0UND=0 independent variables are unbounded 

KB0UND=1 independent variables are non-negative 

KB0UND=2 independent variables are bounded by specific values 

KCLIMB=0 modified objective function optimized via Hooke and Jeeves pattern 
search 
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TABLE 1 (cont'd) 

NOMENCLATURE IN MAIN PROGRAM AND SUBROUTINE OBOECT 


KGRID=0 

KGRID>0 

KK0FF=0 

KK0FF>0 

KMIN=0 

KMIN=1 

KSCALE=0 

KSCALE=1 

KSCALE=2 

KSENS=0 

KSENS>0 


program uses fine-grid hilldimbing convergence criterion at all 
penalty levels 

program uses fine-grid hillclimbing convergence criterion at early 
penalty levels 

automatic termination if Z not uniquely determined by X-vector 
computation continues if Z not uniquely determined by X-vector 
program maximizing the objective function 
program minimizing the objective function 
constraints scaled automatically 
constraints not scaled 

constraints not scaled during a penalty level due to detection of 
suspected constraint ridge by automatic testing procedure 

final values of partial derivatives are not calculated 

partial derivatives of objective function are calculated after 
optimization has been completed 


NCHECK=0 program terminates without a final random check 

NCHECK>0 program generates NCHECK thousand random check points after the 
optimum 


NRAND=0 program performs one search from specified starting point 
NRAND>0 program performs NRAND searches from specified starting point 
Definition of Variables Used in OBJECT 
AREA: area of fuel cell, cm^/cell 

Cll, C12, C13, C14; coefficient used in calculating equation convergence of 

water shift reaction 


CO: current density, A/cm^ 

CNl: Q*A/CMIN in heat exchanger 
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TABLE 1 (cont'd) 

NOMENCLATURE IN MAIN PROGRAM AND SUBROUTINE OBJECT 

CKl, CK2: coefficient used in calculating kinetic equation of demethanation 

EXT: extra air in burner 

NST: number of stacks in fuel cell 

PH20: mole fraction of H 2 O 

P02: mole fraction of O 2 

TAT: ambient temperature, K 

TOPFC: operating temperature of fuel cell, K 

Ul: total heat transfer coefficient of heat exchangers (assume the same 

U value for all exchangers and condenser), caWm^-hr-c 
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The nomenclature of this subroutine is listed in Table 1. The program control 
variables are as follows; 

K*l: read optional input and return 
K*2: evaluate objective function and return 

K»3; evaluate objective function, write optional output and return 
K®4; scale constraints, evaluate objective function, and return 
K»5: write optional output and return 

K=6: evaluate objective function and unsealed SIGMA and return 
ISCALE=1: constraints as yet unsealed 
ISCALE=2: constraints already scaled 

ISCALE=3; constraints scaling in progress 

Subroutine PATTRN 

This subroutine uses the Hooke and Jeeves algorithm (Ref, 4), known as 
pattern search, to do multidimensional search for the better point. This 
technique alternates sequences of local exploratory moves with extrapolations. 
The basis for the method is the intuitive presumption that a strategy which was 
successful in the past will be successful in the future. 

The computation begins by selecting an initial exploratory point Xp and 

evaluating the objective function y(xj^). Then we increase Xj^ to x^^ + and 
evaluate y at this new point. If y shows improvement as a result of this move, 

then we let this new point be the next exploratory point, If y does not 

improve at x^^ + Ep we then try Xj^ - E^. In the event that y does not show 

improvement at both cases, we set equal to the initial exploratory point 



From we proceed to change by + Kg and, thus establish the position 

of the next exploratory point x , as above. The procedure is repeated until 

12 

all the independent variables have been perturbed by some distances. We shall 
refer to this terminal point as a base point, x^. 

A pattern move is now made, 

Xgs x^ + 2 (Xj^ - xj) (19) 

Another sequence of exploratory moves is then carried out. The best value of 
y found during this exploratory process is compared with the previous base 
point. If it is an improvement, it is accepted as new base point. 

The next pattern move is 

Xg = x^ + 2 (Xg - (20) 

and after that each pattern move (except the first) is made in accordance with 
the expression 

^i+1 = ^i-1 - X ._^) (21) 

In case that a pattern move is not successful, then x^-+j^ is set equal to x^ 
so that a new exploratory sequence is initiated from the previous base point. 
The successive explore-and-extrapolate strategy is again applied as before, 
providing that the objective function again continues to improve. 

The flow chart of this subroutine is shown in Figure 4. 
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Subroutine SCALE 


This subroutine computes the multiplicative scaling constants SC(I) such 
that the sums of the magnitudes of the partials of the constraints are roughly 
equal. The constraints in any one problem are often concerned with quantities 
of entirely different magnitudes. In one constraint, one may be dealing with 
values of the magnitude of 10®, and in another with quantities having a 
magnitude of 1. In this form, change in the former constraint will completely 
mask any changes in the latter. When this situation arises, meaningful results 
cannot be obtained unless one "scales" all constraints. This subroutine uses 
the Keefer and Gottfried (Ref. 7) method to scale constraints. 

The flow chart of this subroutine is shown in Figure 5. 

3.2 Program Operation 

The only program input consists of NAMELIST datasets which are read in a 
specified order. The first NAMELIST set is called CON, and contains the 
general information of the optimization procedure and the information of the 
specified problem to be solved. The order of input data inside one NAMELIST 
is unimportant. 

The second set, INIT, contains the numerical values for the flow rate of 
the input fuel, the ambient temperature, the operating conditions of heat 
exchangers, and the operating conditions in the fuel cell stacks. 
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Figure 4 Chart of Subroutine PATIEN 
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Figure 4 continued 





The third set, CONI, carries the constants used in evaluating the reaction 
rate of demethanation reaction and the equilibrium constant of the water shift 
reaction in the fuel processor. 

The fourth set, C0N2, contains the coefficients used in estimating the 
cost of the power inverter. 

The last NAMELIST dataset, SI, gives the initial guesses for the flow rates 
and temperatures of specified flowsheet streams {Figure 1). In addition, the 
bounds and initial guesses of decision variables are also listed in SI. 

More detailed descriptions of these input data are shown in Table 2, along 
with their dimensional units, and the numerical values used in the sample run. 

Besides the input data, the values of some variables have been fixed in 
the program as the initial assumption, which is listed in Table 3. 

3.3 Sample Run 

The object of the sample problem is to calculate the optimal values of 
hydrogen utilization (x^) in the fuel cell stack, the steam to carbon ratio 
(S/C, X 2 ) in the reformer, and the number of cell plates per stack (x^) in 
order to minimize equipment capital cost of the 100 kW (AC) PAFC powerplant 
system. The objective function and constraints are listed in Chapter li. The 
input data for this sample problem is shown in Table 2. To avoid carbon 
formation inside the reformer, the minimal value of X 2 was assigned a higher 
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TABLE 2 


INPUT DATA FOR OPTIMIZATION COMPUTER MODEL (SAMPLE RUN) 


N; 

number of independent variables 

3 

OEQ: 

number of functional equality constraints 

0 

OINEQ: 

number of functional inequality constraints 

10 

IPMAX; 

maximum number of penalty levels 

25 

ALPHA: 

scaling factor between SIGMA2 and SIGMAl 

1, 

ARAND: 

initial parameter for the random number generator 

0.314159167 

CFAC: 

proportionally constant used in calculating PNLTY 

0.01 

CONST: 

penalty multiplier 

8. 

CONVG: 

penalty level convergence criterion 

0.0001 

DELTA: 

finite-difference increment 

1.0-10 

DX: 

length of a gradient move 

10. 

EFAC: 

search increment reduction factor 

0.001 

EPSLN: 

gradient vector convergence factor 

l.D-05 

GRADX: 

gradient vector convergence criterion 

0. 

THETA: 

maximum time allowed to evaluate the objective func- 
tion when calculating the gradient vector by central 
differences 

0.1 

I0UT=0 

output follows final optimization 


KB0UND=2 

independent variables are bounded by specific values 


KCLIMB=0 

modified objective function optimized via Hooke and Jeeves pattern 
search 

KGRID=1 

program uses fine-grid hillclimbing convergence criterion 
penalty levels 

at early 

KK0FF=0 

automatic termination if Z not uniquely determined by X-vector 


program minimizing the objective function 


KMIN=1 


TABLE 2 (cont'd) 

INPUT DATA FOR OPTIMIZATION COMPUTER MODEL (SAMPLE RUN) 


KSCALE=0 constraints scaled automatically 

KSENS=0 final values of partial derivatives are not calculated 

NCHECK=0 program terminates without a final random check 

NRAND=0 program performs one search from specified starting point 

AREA; area of fuel cell, cm^/cell 1316.13 

Gll, C12, C13, C14; coeffiecients used in calculating equation convergence of 

water shift reaction; -9.0283D+8, 3.5603D+6, 4.3662D+3, 
~3.0526 

p 

CD: current density, A/cm 0.1 

CNl: U*A/CMIN in heat exchanger 1.3 

CKl, CK2: coefficients used in calculating kinetic equation of demathanation 
reaction; 127., -4417.7 

DDl: amount of input methane, g-mole/hour 1297.3 

EXT: exta air in burner, % 100. 

NST: number of stacks in fuel cell 6 

PH 2 O: mole fraction of H 2 O 0.095 

PO 2 : mole fraction of O 2 0.1413 

TAT: ambient temperature, K 298. 

TOPFC; operating temperature of fuel cell, K 463. 

Ul: heat transfer coefficient of heat exchangers, 

cal./m2-hr-K 48825.1016 

CCl, CC2: power inverter cost data 

CCl 50., 164., 203., 248., 304., 366., 433., 528., 657., 920., 1560., 

2810., 1000000000 

CC2 200., 160., 150., 140., 130., 120., 110., 100., 90., 80., 70., 60., 

50 
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TABLE 2 (cont'd) 

INPUT DATA FOR OPTIMIZATION COMPUTER MODEL (SAMPLE RUN) 






min. 

max. 

i n i t . 

X(l): 

utilization of H 2 in 

the fuel cell 

0.7 

0.95 

0.89875 

X(2): 

steam/ methane ratio 



1.7 

3. 

1.7192 

X(3): 

number of cell plates per stack 

200. 

250. 

225.38 

S(l): 

flow rate of CH^ in 

steam 7; Ib-mole/hr 



0.51439 

S(2): 

CO 

II 




1.7525 

S(10): 

II 


9 



1.3551 

S(7): 

II 


13 



0.018455 

S(5): 

COj 


7 



0.63024 

S(17): 

M 


9 



0.98925 

S(9): 

II 


13 



2.3257 

S(8): 

H 2 


7 



7.6278 

5(12): 

II 


9 



8.0222 

5(3): 

M 


13 



9.3582 

5(11): 

H 20 


7 



2.5643 

5(13): 

M 


9 



2.1926 

5(4): 

II 


13 



0.8295 

5(14): 

temperature of stream 8; 

K 



879.31 

5(16): 

II 
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637.04 

5(18): 

M 

5 




1096.6 

5(16): 

II 

18 




1197.7 

5(6): 

It 

16 




772.8 

5(20): 

amount of catalyst in reformer section; lb 



60. 
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TABLE 3 


ASSUMED PARAMETERS 

1. input of fuel (methane) - 1297.3 g-mole/hour 

2. designed AC output - 100 kW 

3. average current density on the fuel cell plate - 0.1 A/cm^ 

4. operating pressure in the reformer and fuel cell stack - 1 atm. 

5. utilization of O 2 in the fuel cell stack - 50 

6. catalyst factor (--utilization x surface area x loading) - 45 

7. dimension of cell plate - 17x12 in^ 

8. ADT (approach differential temperature) of shift rxn. in the reformer 
and shift converters - 15 K 

9. utilization of air in the burner - 50 

10. the bulk density of catalyst in the reformer - 30 1b/ft^ volume of 
reformer 

11. temperature of inlet steam - 391 K 

12. operating temperature of fuel cell stack -463K 
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value than the minimum steam to carbon ratio (Refs. 8 and 9). There are no 
further restrictions for one to assign the maximum and minimum of variables, 
but the reasonable values will result in less computation time. The values 
chosen in this sample problem are listed in Table 2 under NAMELIST dataset 
SI. Some of the operating conditions of the powerplant model used as input to 
this optimization procedure are derived from the steady state performance code 
(Ref. 10). 

The output of this sample problem is shown in Figure 7. First, the 
program reprints all of the input data. The optimal solutions of three design 
variables are given under the policy vector, and the minimal cost is the value 
of Y. The sample output are given below: 

Hydrogen utilization (xj^) 90.6 

Steam to carbon ratio (X2) 1.7 

Number of cell plates per stack (X3) 226 

Minimal equipment capital cost of 100 kW (AC) PAFC powerplant $33534 

The required CPU time to run this sample problem is less than 1 minute on 
IBM/370. 

3.4 Discussion 

It can be noted that the optimal values are some differnt from the usual 
design. The optimal hydrogen utilization in the stack is higher than the 
designed one (0.906 and 0.75-0.85, respectively) , and the optimal S/C (1.7) 
is close to the minimum value. 
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ORIGINAL IS 
OF POOR QUALITV 


ICON JINEQ=10» 

JEQ=0, * 

KB0UND=2, 

KMIH=1, 

NRAND=0, 

KCLIMB=0, 

KKOFF=0» 

MPEN=0, 

IOUT=0, 

NCHECK=0, 

IGRAD=1, 

K5ENS=0, 

KGRID=1, 

KSCALE=0, 

IPMAX=25, 

ALPHA=l.D+00, 

ARAND-0.314159267» 

BOUNDS6.023D+23, 

CFAC=l.D-02, 

PNLTY=l.D+00, 

CONST=8.D+00, 

DELTA=1.D-10, 

DX=l.D+02, 

EFAC=l.D-03, 

EPSLN=l.D-05, 

THETA=1.D-01, 

TIME=0. , 

C0NVG=1 .0-4, 

GRADX=0. , 

lEND 

UNIT DD1=1297.3,TAT=298.,P02=0.1413,PH20=0.095,TOPFC=463. ,CN1=1.3, 

U1 =48825 , 1016 , EXT = 100 . ,NST=6 ,CD=0 .100, AREA =1316. 13, T26'<"39i . , 

SEND 

SCGNl CK1=127. ,CK2=“44 17. 7, Cll=-9. 02830+8,012=3. 56 03D+6,C13=4. 36620+3, 
C14=-3.0526, 

SEND 

IC0N2 CC1=50., 164., 203., 248., 304., 366., 433., 528., 657., 920., 1560., 

2810. , 1000000000. , 

CC2=200,, 160., 150., 140., 130., 120., 110., 100., 90., 80., 70., 60., 

50. , 

SEND 

SSI XMIN(1)= 0.7, 

XMIN(2)= 1.7, 

XriIN(3)= 200.;. 

XMAXC1)= 0.95, 

XMAX(2)= 3.0, 

XMAX(3)= 250., 

X(1)=0. 89875, 

X(2)=1.7192, 

X(3)=225.38, 

S(1)=0. 51439, 

5C2)=1.7525, 

S(3)=9.3582, 

S(4)=0.8295, 

5C5)=Q. 63024, 

S(6)=772.8, 

5(7)=0. 018455, 

S(8)=7.6278, 

S(9)=2.3257, 

S(10)=1.3551, 

S(ll)=2.5643, 

S(12)=8.0222, 

S(13)=2.1926, 

S(14)=879.31, 

S(15)=1197.7, 

S(16)=637.04, 

S(17)=0.98925, 

S(18)=1096.6, 

S(20)=60. , 

SEND 


Figure 6 Sample Input Data 
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Figure 7 Sample Computer Run 
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Figure 7 (cont'd) Sample Computer Run 


This program developed herein can help set the values of operation 
parameters in order to optimize the cost of a given plant. As an example, it 
might be used to perform the analysis to determine the tradeoff between 
operating a powerplant at a higher pressure and temperature, which improves 
power section efficiency, but also accelerates power section replacement 
schedule, and operating at lower temperature and pressure. These design and 
operating options are very important in establishing technology goals for the 
phosphoric acid fuel cell system. 

» 

3.5 Further Developments 

There are two other optimization programs embedded in the overall strategy, 
which is not included in this manual; refer to Refs. 8 and 9. One subsystem 
program determines the minimum steam to carbon ratio for the reformer. This 
lowest permitted ratio value is a critical constraint for the total optimiza- 
tion model, not only as an operational danger point to the reformer catalyst, 
but also as cost factor to be minimized in the reaction water recovery heat 
exchanger subsystem. The other subsystem optimization program analyzes the 
thermal management heat exchanger network, and can improve that subsystem prior 
to the main optimization analysis of the total system. 
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LISTING OF THE OPTIMIZATION COMPUTER MODEL 
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